Correlated mesoscopic fluctuations in integer quantum Hall transitions 
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We investigate the origin of the resistance fluctuations of mesoscopic samples, near transitions 
between Quantum Hall plateaus. These fluctuations have been recently observed experimentally by 
E. Peled et al. [Phys. Rev. Lett. 90, 246802 (2003); ibid 90, 236802 (2003); Phys. Rev. B 69, 
241305(R) (2004)]. We perform realistic first-principles simulations using a six-terminal geometry 
and sample sizes similar to those of real devices, to model the actual experiment. We present the 
theory and implementation of these simulations, which are based on the linear response theory 
for non-interacting electrons. The Hall and longitudinal resistances extracted from the Landaucr 
formula exhibit all the observed experimental features. We give a unified explanation for the three 
regimes with distinct types of fluctuations observed experimentally, based on a simple generalization 
of the Landauer-Biittiker model. The transport is shown to be determined by the interplay between 
tunneling and chiral currents. We identify the central part of the transition, at intermediate filling 
factors, as the critical region where the localization length is larger than the sample size. 

PACS numbers: 73.43.-f, 73.23.-b, 71.30.+h 



I. INTRODUCTION 

The integer quantum Hall effect 1 (IQHE) is one 
of very few instances where quantum effects are dra- 
matically manifested in the macroscopic world: if a 
two-dimensional electron system (2DES) is placed in 
a large perpendicular magnetic field B, its Hall re- 
sistance is quantized with high accuracy to h/(ne 2 ), 
where n > is an integer. This quantization has been 
explained, within the non-interacting electron approxi- 
mation, as a bulk effect^ as a topological invariant? 3 , 
and as an edge effect^ The plateau-to-plateau tran- 
sitions of the Hall resistance, which are accompanied 
by peaks in the longitudinal resistance, are understood 
as localization-delocalization transitions with a universal 
scaling relationiSiLSi2ii£ 

A different route to observing quantum effects in con- 
densed matter physics, is to reduce the size of the device 
to the so-called mesoscopic regime. If the size of the de- 
vice is smaller than the phase coherence length, quantum 
interference between different paths leads to effects such 
as random universal conductance fluctuations! 11 ! 12 ' 13 

What happens if the two are combined? The answer 
was provided by recent IQHE experiments performed on 
a mesoscopic sample ^iiSiiSi The expected quantum me- 
chanical interference of electronic wave- functions indeed 
causes reproducible patterns of fluctuations in the resis- 
tances near plateau-to-plateau transitions, as the mag- 
netic field is varied (in macroscopic samples all resis- 
tances vary smoothly with B). Interestingly, the fluc- 
tuations of resistances measured with different combina- 
tions of electrical contacts are not independent; instead, 
various non-trivial correlations were observed. 

In this paper, we investigate the mesoscopic IQHE 
within the non-interacting electron approximation. Our 
numerical simulations provide the basis for a unified ex- 
planation of these novel experimental observations, based 



on existing theories of the IQHE and mesoscopic trans- 
port (some of these results have been summarized in 
Ref. We begin the paper with a brief review of 

the experimental observations, in Sec.[H] In Sec. IIIII we 
describe in detail the model we use, and our approach 
to solve the scattering problem that allows us to com- 
pute the conductance matrix from the multi-terminal 
Landauer formula. Our numerical results, which exhibit 
all the experimentally observed symmetries and correla- 
tions, are presented in Sec. II VI In Sec.[V]we explain how 
these various correlations and symmetries arise. We an- 
alyze the allowed structure of the conductance matrix, 
and show how it directly determines the possible correla- 
tions between fluctuations of various resistances. Finally, 
Sec . I VII contains the summary and a discussion of the rel- 
evance of these results. 



II. SUMMARY OF EXPERIMENTAL RESULTS 

The experimental sample is an InGaAs/InAlAs het- 
erostructure with impurities (Indium atoms) distributed 
throughout. As a result, the mobility is rather low (frac- 
tional quantum Hall effect is not observed) . Even for the 
IQHE, only the first few plateau-to-plateau transitions 
can be identified. The experimental setup is sketched 
in Fig. ^ The Hall bar is connected to six terminals, 
indexed 1 to 6. The sample is considered mesoscopic be- 
cause the characteristic size of its central region between 
the 4 central terminals, of about 2^tm, is comparable to 
the phase coherent length w 2 — 3/LtmiMiiSii& 

As customary, we denote the various resistances as 



n /Iij, where V mn = V n - V m 



is the voltage 

difference measured between terminals n and m, when 
the current is injected in the sample through terminal i 
and extracted through terminal j. For example, to mea- 
sure the Hall and longitudinal resistances, a small current 
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W=2-4nm 



FIG. 1: (color online) Typical six-terminal Hall bar used in 
experiments. In this configuration, the Hall resistance Ry^m 
and the longitudinal resistance ^14,23 are measured. 



D: Under the reversal of the magnetic field B — » —B, 
the longitudinal resistances verify the symmetrytiS 



#14.23 (#) 



#14,65 (~ B )< 



(3) 



to high accuracy, although Rx^zi^) ^ #U23( — #) an d 
i?f 4 65 (B) ^ i?f 465 (— S) except in regime A3, on the 
high-zy side of the transition, where the two longitudinal 
resistances have identical fluctuation patterns. 

In the following, we show that all these observations 
can be explained within the non-interacting electron ap- 
proximation, using a combination of ideas about the 
IQHE and mesoscopic transport. 



then measure two longitudinal resistances i?f 4 23, *H4 65; 
and two Hall resistances i?i4. 6 2; #14,53 



/ is fed into terminal 1 and collected from terminal 4; no 
net currents are flowing in the other 4 terminals. One can 

R{ 

For any macro- 
scopic sample, the two resistances of each pair are equal. 
For a mesoscopic sample, however, each resistance ex- 
hibits a different fluctuation pattern>iii^i& 

Another quantity of interest is the "two-terminal re- 
sistance" i?2t = #63,63- Experimentally, it is measured 
by sending the current from terminal 6 to terminal 3, 
and measuring the voltage difference between terminal 6 
and terminal 3. Because of the two-point measurement, 
a contact resistance is subtracted from the raw data^ 
Summary of experimental findings: 
A: the IQHE transition from the n th to the n + 1 st 
plateau, Vn > 0, has three distinct regimes»^*i& 

Al: on the high-i? (low filling factor v) side of the 
transition, the Hall and the longitudinal resistances have 
distinct fluctuation patterns, which however are corre- 
lated such that: 



#14,62 + #14,23 



R 



H 

14,53 



R 



h 



14,65 



(1) 



A2: in the central part of the transition, the Hall and 
the longitudinal resistances still exhibit distinct fluctua- 
tion patterns, but Eq. is no longer satisfied. 

A3: on the low-fi- (high filling factor v) side of the tran- 
sition, the Hall resistances become quantized to the ex- 
pected plateau value R^ 62 = #{4 53 = h/[(n+l)e 2 ]. The 
two longitudinal resistances exhibit significant, but this 
time identical fluctuation patterns, with a maximum am- 
plitude h/[n(n + l)e 2 ], i.e. the difference between quan- 
tized Hall resistances of the two neighboring plateaus. 

B: the transition inside the lowest Landau level n = 
shows the equivalent of A3: at high-z/ one finds the regime 
where i?i4. 62 = #14,53 = h/e 2 , while the two longitudi- 
nal resistances have identical fluctuation patterns with 
extremely large amplitudes. 14 Regimes Al and A2 are 
replaced here by the transition to the insulating phase, 
where all four resistances increase rapidly with decreas- 
ing filling factor. 

C: throughout each IQHE transition, the following 
identity is found to hold to high accuracy*^ 



#14 fi9 + #14 



R 



H 

14,53 



R 



14,65 



R 



63,63- 



(2) 



III. NUMERICAL SIMULATIONS 
A. The conductance matrix and the resistances 

The linear response function of the six-terminal Hall 
bar is the 6x6 conductance matrix g that defines the re- 
lationship between the currents flowing through the var- 
ious leads, and their voltages: I a = J2p QapVp. Here, I a 
is the current in the lead a — 1, . . . , 6. We use the con- 
vention / > (/ < 0) for currents coming out of (flowing 
into) the sample. V a is the voltage of the lead a. 

Knowledge of the conductance matrix allows us to cal- 
culate the various resistances in a straightforward way. In 
the conventional setup, the current flows from lead 1 to 
lead 4 and therefore J 14 = (-J I 0) T . Without 
loss of generality, we set I — 1 (this is a linear response 
theory) and V4 = (the voltage differences are not af- 
fected if one of the terminals is grounded) . After solving 
the 6x6 equation gV = /14 for the other 5 voltages, we 
can directly compute the various resistances from their 
definitions: 



#14,23 — Vz 
#14,62 =Ve — V2, 



V3, #14.65 



R 



H 

14,53 



= u 6 - 
= v 5 



v 5 , 
-v 3 . 



(4a) 
(4b) 



The two-terminal resistance is determined similarly. In 
this case, the electric currents flow from lead 6 to lead 
3: J 63 = (0 I ~I) T . After solving the equa- 
tion gV = Iq3 for the corresponding voltage distribution, 
we can calculate immediately the two-terminal resistance 
(assuming again that 1=1): 



i?2t — # 



63,63 



= V*-V* 



(5) 



Clearly, any other measurement can be simulated in a 
similar fashion within this formalism. All resistances are 
functions of various elements of the conductance matrix. 

At T = 0, the off-diagonal elements of the conductance 
matrix g are calculated^ by solving a multi-channel scat- 
tering problem, according to the Landauer-Buttiker for- 
mula: 



9a,l3^a{E F ) = — y \t a i,f3j{E F )Y 



(6) 
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Here, t a i i pj(Ep) is the amplitude of probability (trans- 
mission amplitude) that an electron with the Fermi en- 
ergy Ep, which is injected into the sample through the 
j th transverse channel of lead /3, will scatter out of the 
sample into the i th transverse channel of lead a. The 
sum over all the transverse channels of the two leads sim- 
ply gives the total probability pp^ a (Ep) for the electron 
with Ep to scatter from lead /3 into lead a, and thus: 



(7) 



The diagonal elements of the conductance matrix can 
then be calculated using the restrictions imposed by 
charge conservation and gauge invariance: 18 ^ Q g a p = 
J2a9a/3 — 0- The diagonal elements are thus 



9o 



(8) 



The finite-temperature generalization in the linear 
regime is straightforward: 



g a ,p(ji,T) 



dEg ai0 (E) 



df(E) 



dE 



(9) 



where f(E) — [exp(/3(E — /i)) + l] 1 is the Fermi distribu- 
tion, where (5 — l/(fcsT) and /i is the chemical potential. 

The strategy for the numerical simulations of the IQHE 
in mesoscopic samples is thus apparent. We have to solve 
a complicated, multi-channel scattering problem to find 
all the amplitudes of transmission through the Hall bar. 
This allows us to calculate the conductance matrix and 
therefore the various resistances. This is repeated for 
many values of the Fermi energy (corresponding to dif- 
ferent filling factors v of the Landau level) to obtain the 
the traces of the resistances as v changes. We now de- 
scribe the details of our simulations. 



B. The model 

In this section we describe the model we use to nu- 
merically simulate IQHE transitions in mesoscopic sam- 
ples. This comprises the sample, the confining potential, 
the disorder potential, the external leads and their con- 
tacts to the sample. As already stated, we assume non- 
interacting electrons - the usual approximation when 
treating the IQHE. For numerical convenience, we ignore 
Landau level (LL) mixing, although the whole formalism 
can be straightforwardly extended to include it. LL mix- 
ing can be ignored for large magnetic fields (the case of 
interest to us), when the cyclotron energy is the largest 
energy scale in the problem. In the following, we take 
the effective mass of the electron in the 2DES to be that 
of GaAs/AlAs heterostructures, m* — 0.067m e , and use 
— e, e > for the electron charge. 



1. The sample 

In choosing the sample geometry, we face an appar- 
ent problem. On one hand, when working with Landau 
levels, it is convenient to consider a rectangular sample 
with periodic boundary conditions along one axis, since 
in this case the eigenstates of each LL are known analyti- 
cal functions. This simplifies considerably the calculation 
of various matrix elements. On the other hand, however, 
a Hall bar with different contacts on all sides, such as 
sketched in Fig. ^ is obviously not consistent with peri- 
odic boundary conditions. 

Our solution is to start with an area larger than the 
Hall bar itself. For this enlarged area we impose cyclic 
boundary conditions to generate a convenient basis for 
the Hilbert subspace of each LL. We then add a confining 
potential to "carve out" the Hall bar from this larger 
area. This potential is described in the next subsection. 

Consider, then, 2D electrons in an area of size L x x L yi 
with periodic boundary conditions in the y-direction, 
and placed in a perpendicular magnetic field B. The 
spectrum consists of Landau levels of energy E na — 
huj c (n + 1/2) — guBBa, where n = 0,1,2,... is the LL 
index, and a = ±1/2 is the z-axis spin projection. The 
cyclotron energy is hu> c = eB/m*. 

If we use the Landau gauge A = (0, Bx, 0), the eigen- 
functions for the (n, a) Landau level are>i£ 



(rln, X,a) = 



y/2 n n\y/wl' 



(10) 



where I = y ^ is the magnetic length, H n (x) are Her- 

mite polynomials and Xa- are the z-spin eigenstates. The 
cyclic boundary condition in the y-direction require that 
Xj = j2irl 2 /L y , j E Z. Xj, the guiding center, charac- 
terizes the location at which individual basis states are 
centered on the x-axis [see Eq. <|10[) ]. Since —L x /2 < 
Xj < L x /2, it follows that -N/2 <j< N/2, and the de- 
generacy of each spin-polarized LL is N = L x L y /(2iil 2 ). 
As we show later, TV roughly defines the size of the ma- 
trix involved in solving the scattering problem (several 
more states on the contacts have to be included as well) . 

In the simulations shown here, we use L x = L y = 4/im. 
For a field of several Tesla, this leads to a degeneracy N ~ 
10 4 , which can be easily handled by a generic PC cluster 
in a reasonable amount of time. The main problem is 
that our sample is shorter than the experimental sample 
(4/Ltm instead of ^20/^m). This means that terminals 1 
and 4 (see Fig. ^| are much closer to the other 4 terminals 
than in reality. We fix this problem with the aid of the 
confining potential, as described in the next subsection. 



2. The confining potential 

As stated before, we add a confining potential V c {r) to 
define the Hall bar from the larger area L x x L y spanned 



FIG. 2: The confining potential V c (r) that we use. It varies 
smoothly near the y-edges, and has triangular barriers in the 
corners. See text for details. 

by the LL Hilbert subspace. We have tested several 
functional forms, to see which are physically reasonable. 
Some of the choices considered and their pros and cons 
are discussed in Appendix D 

The confining potential used in simulations is shown in 
Fig. [3 It is an odd function V c (x,y) = —V c {x,—y) and 
such that V c (x,—L y /2) = 0. The Hall bar is the region 
—Ly/2 < y < 0/2 and —L x /2 < x < L x /2, where the 
confining potential is: 

(1 _ e »A)d - e — ^r^) 

KM = "W !_ e -W" X 

J] 9(\x- vL x /2\,\y + nL y /2\), (11) 

v = ±l 
f = 0,1 

The first line describes the main features: the confin- 
ing potential is approximatively equal to — Vg ap inside 
the Hall bar. We use V gap = 10 meV in our simula- 
tions, which is much larger than the amplitude of the 
disorder potential (see below). As a result, for Fermi 
energies Ep < 0, the electrons are confined to the re- 
gion [—L x /2,L x /2] x [— Ly/2,0] which defines the Hall 
bar to which contacts will be attached. Near the Hall 
bar y-edges at y — 0, —L y /2, the confining potential rises 
smoothly to zero, on a length-scale A (chosen to be 40 nm 
in our simulations; see discussion in Appendix^). The sec- 
ond line describes the triangular potential barriers added 
in the four corners of the Hall bar: 

g{x,y) = sm 2 [(x/l x + y/l y )ir/2]Q(l - x/l x - y/l y ) + 
Q{x/l x + y/l v - 1), (12) 

where l x and l y — L y /6 define their spatial extension. 
These barriers help isolate the contacts 1 and 4 from the 
other contacts, and thus compensate for the shortness of 
our sample. The alternative is to use a longer sample; 
this, however, requires significantly more CPU time. 

The symmetric choice we made for the confining po- 
tential is not necessary; however, it is advantageous for 
several reasons. If we take the disorder potential such 



FIG. 3: (color online) Total confining plus disorder potential 
inside the Hall bar, —Ly/2 < y < (compare with Fig. El - 

that Vd(x, — y) = —Vd(x,y) as well, then the particle- 
hole symmetry of the total potential for the whole sam- 
ple means that all quantities must be symmetric with 
respect to Ep = 0. This allows us to check and calibrate 
the numerical procedure, as discussed in Appendix D 

Another advantage of a symmetric choice is that it 
allows us to easily define the Hall bar filling factor for 
Ep < 0, as being twice the filling factor of the LL (same 
number of filled states, but half the area). Strictly speak- 
ing, this is an underestimate. At the Fermi energies 
Ep ~ — Vg a p of interest to us (see below), the effective 
area of the Hall bar is somewhat smaller, because of the 
shape of the confining potential. 



3. The disorder potential 

The experimental sample contains many In impurities 
which serve as short range scattering centers. We sim- 
ulate a short-range disorder potential Vdif) by adding 
18000 Gaussian scatterers randomly placed inside the 
Hall bar area [—2, 2]/xmX [—2, 0]/im. Each Gaussian scat- 
tering potential is of the form Ae~ r / d , where d is uni- 
formly distributed in the interval [0.01, 0.03) /im, and A 
is uniformly distributed between [—0.6, 0.6] meV. 

A typical disorder realization within the Hall bar is 
shown in Fig. |3 added to the confining potential for the 
Hall bar side of the total sample. The disorder is very 
rough and short-range, almost white-noise-like. More- 
over, its amplitude is much smaller than that of the con- 
fining potential, as required in order to have a well de- 
fined Hall bar. Our simulations scan a dense energy grid, 
typically in the range [—10.3, —9.7] meV, close to a filling 
factor v ~ 1/2 for the Hall bar. 



4- The leads and the contacts 

The IQHE is independent of the type of current- 
carrying external leads used and of their contacts to the 
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sample, provided that (i) the leads are reasonably good 
metals and (ii) the contacts allow easy transfer of the 
electrons into and out of the Hall bar. As a result, no de- 
tailed realistic modeling of the leads and of the contacts 
is required. Instead, we can use simple, idealized models 
that satisfy requirements (i) and (ii). 

As in Ref. [lj| we model each external lead as a col- 
lection of independent, semi-infinite, perfectly metallic 
one-dimensional tight-binding chains. Each such chain 
represents a transverse channel of that lead, which car- 
ries currents independently of the other channels in the 
lead. The choice for the hopping matrix t and on-site en- 
ergy eo that define the Hamiltonian of each such chain, 
has been discussed in detail in Ref. Briefly, t must be 
chosen so as to minimize the contact resistance. We use 
t = 0.1 meV, comparable to the major matrix elements 
of the Hamiltonian matrix elements inside the sample. 
To insure that the leads can always carry currents to and 
from the Hall bar, we use a "floating" spectrum, i.e. we 
set eo = Ep that we investigate. With these choices, 
requirement (i) is trivially satisfied. 

The contacts are represented by matrix elements con- 
necting states on the leads (tight-binding sites) to contact 
states in the Hall bar, which are appropriate linear com- 
binations of the LL basis states. The simplest choice is to 
add to the total Hamiltonian a hopping term of the form 
— t(c ! Q d + h.c.) for each chain (transverse channel). Here, 
Cq is the creation operator for the last site of the chain, 
while d is the annihilation operator for its corresponding 
contact state. 

The contact states are chosen so as to satisfy the fol- 
lowing general requirements: 

(1) Each contact state must be localized close to the 
region of the Hall bar where the corresponding lead is 
connected. 

(2) Different leads have orthogonal contact states, but 
channels of the same lead may have non-orthogonal con- 
tact states. This is because the direct connection (short 
circuit) between different leads is forbidden. 

(3) For numerical efficiency, it is convenient to preserve 
the sparsity of the total Hamiltonian. We therefore re- 
quire that a contact state is a linear combination of a 
small subset of the N LL basis states. 

For leads 1 and 4 (the current source and drain, see 
Fig. HJ, the contact states must be localized close to the 
x = ±L x /2 ends of the sample. The LL basis states 
\n,Xj, a) with j ~ ±N/2 are exactly such states. For 
simplicity, we take each of them to be a contact state, 
as in Ref. 0. Specifically, we assume that leads 1 and 
4 have N c transverse channels each (we use N c — 100 in 
the simulations shown here) . The m th transverse channel 
of the lead 1 (4) is coupled through a hopping term to 
the LL basis state \n,Xj,a) with j = —N/2 + m, and 
respectively j = N/2 — m, (V) 1 < m < N c . With this 
choice, all contact states are orthogonal to one another, 
and all three requirements listed above are satisfied. This 
means that electrons can be injected/removed from the 
Hall bar anywhere within a strip of width N c 2irl 2 / L y 




FIG. 4: (color online) Schematic drawing of the probability 
density of a contact state on the y-edge. It overlaps with the 
entire region where the confining potential rises from — Vg ap 
to Igap, so that it can inject /collect electrons of any Fermi 
energy from the Hall bar. 

from the x = ±L x /2 edges, through contacts 1 and 4. 

For the voltage probes (leads 2, 3, 5 and 6) we need 
a different strategy, since they are located on the y = 
—Ly/2, respectively y — edges of the Hall bar. The 
contact states are chosen to be of the form: 

(^vacuum) = A^ e~ 'V +im6 \n,X m ,a) , (13) 

m 

where a is a parameter that controls the x and y extension 
of the wave-function, 9 = 0, tt for y = 0,—L y /2 respec- 
tively, and A is the normalization factor. As a changes, 
the shape of this wavefunction changes from a strip in 
the s-direction to a strip in the y-direction with fixed 
total area. We choose a — 0.21 so that the probability 
distribution of the contact states overlaps with the region 
of smooth variation of the confining potential near the y- 
edges. In this case, the profile of the contact states is an 
ellipse with its longer axis in the y-direction (see Fig.^J. 

In Eq. I|13fl , x c specifies the location of the contact state 
along the y-edge. We use Ac/4 (25, in these simulations) 
tight-binding chains for each of these leads, and corre- 
spondingly N c /A different contact states. For leads 3 and 
5, the centers x c of their contact states are distributed 
uniformly in the range [1, 1.5]^m, while for leads 2 and 
6 the range is [—1.5, — l]/xm. The summation in Eq. i|13|) 
is truncated so that no pair of neighboring contact states 
contain the same LL basis states in their expansions. The 
contact states on the y — —L y /2 edge (for leads 2 and 3), 
are magnetic translations of the corresponding states of 
leads 5 and 6. The alternating phase factor e lm7r insures 
that contact states on different y-edges are orthogonal to 
one another, so that short-circuits are avoided. 

C. Solving the scattering problem 

The total Hamiltonian is the sum of the sample Hamil- 
tonian, the Hamiltonians for the 6 leads, and the terms 
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describing the contacts between sample and leads: 



it — ^sample 



6 

rl Iead 

a=l 



He 



(14) 



The numerical results presented here are obtained in the 
lowest Landau level (LLL) with n = (the spin projec- 
tion is irrelevant. Calculations in higher LL proceed in a 
similar fashion). For simplicity, in the following we de- 
note the LLL basis states by \Xj) - 

ac 

LLL, the Hamiltonian of the sample is: 



cj I vacuum; = \n = 
0,Xj,a). In the Hilbert subspace of the spin-polarized 



N/2 N/2 

H sam ple= E (XilVc+VdlXj)-^, (15) 

i=-N/2 j=—N/2 

where N = L x L y / (2irl 2 ) is the number of states in a 
LL and an overall constant LLL energy shift Eq has 
been ignored. The matrix elements (X^V^ + Vd\Xj) of 
the confining and disorder potentials are computed as 
described in Ref. Briefly, the idea is that matrix 
elements of a plane- wave are simple analytical functions: 



(X, 



iqr 



x 4 ) = s 



Xi.Xi 



-'-hi 



,2 e 



where Q = \l 2 {q^ + Qy)- (The generalization for higher 
LLs and/or LL mixing is straightforward, see for in- 
stance Eq. (7) in Ref. The strategy, then, is to 
perform a fast Fourier transform of the potentials on a 
grid dense enough to reproduce them with high accu- 
racy while maintaining the sparseness of the Hamilto- 
nian. The matrix element of each Fourier component is 
given by the previous formula, allowing one to efficiently 
compute the matrix elements of the confining and dis- 
order potentials. Note that for higher LLs, the matrix 
elements differ only by a Laguerre polynomial. This is 
not likely to influence the physics in a significant way and 
suggests that numerical results for higher LLs should be 
qualitatively similar to those we present here for the LLL. 

The leads are collections of N" semi-infinite tight- 
binding chains, each describing a transverse channel: 



'head — 



E E h 



(<*)U«) 



a)t c (Q) + hc 

c i,m+l ' ll - c - 



(16) 

Thus, cf^ creates an electron on the site to > of the 
i th transverse channel of lead a. We use the convention 
that site m — is always the end site of the chain. The 
solution we present can be trivially generalized to allow 
for different eo an d t parameters for each channel, as well 
as longer-range hopping. However, since the physics in 
the sample is independent of the lead details (as long as 
the leads can carry currents to and from the sample) the 
simple choice we make should be sufficient. The choices 
for i, €q and number of channels N" for each lead have 
been discussed previously. 



Finally, the contacts are described by matrix elements 
between sites to = of the various channels and their 
corresponding contact states: 



^EERo 

a— 1 i—l 



h.c. 



(17) 



where, as described in the previous section, each contact 
state is a linear combination of LLL basis states, which is 
localized in the appropriate region of the sample. Again, 
generalizations for more complex contacts can be easily 
incorporated, but should not be necessary. 

We now have to solve a scattering problem for an 
electron injected with energy E in the channel j of 
lead /?, which means that we must find an eigenstate 
HlVpj) = E\$ Pj ) of the form 



I* 



0il 



N° 



EEE< 

a— 1 i—l rn—0 



4 («U«)t 

i.m i.m 



N/2 

E < 

j=-N/2 



>3 C 3 



where - with our convention for indexing channel sites 



ihm 



0& 



-tkm 



In other words, there is 



an in-coming wave with unit amplitude in channel (f3j), 
and out-going waves with various transmission ampli- 
tudes in all other channels. The momentum k > is 
such that E = eg — 2tcosk 7 i.e. it is the momentum 
of an electron with energy E propagating on any of the 
tight-binding chains. 

We solve this scattering problem by recasting it as 
a finite-size, inhomogeneous system of linear equations, 
which can be easily solved numerically. The unknowns 
are the amplitudes (f>j to find the electron in various LL 

states in the sample, plus 5 terms (pf^, to = 0, . . . , 4 for 
each channel (clearly, if one knows the wave-function at 
the first few sites of a tight-binding chain, one knows the 
wave-function along the entire chain). This solution is 
explained in detail in Ref. 0, for two leads with multi- 
ple channels. The redistribution of the channels amongst 
more leads is trivial, since it only involves changing some 
of the matrix elements to the appropriate contact states. 

The solution of this inhomogeneous system of linear 
equations, then, directly gives us the transmission coeffi- 
cients t a i t f3j (E) from which we compute the conductance 
matrix [see Eq. (|SJ|]. This is repeated for many Fermi 
energies Ep corresponding to various filling factors v, to 
obtain the dependence g(v) of the various elements of the 
conductance matrix on the filling factor. 



IV. NUMERICAL RESULTS 

A. The conductance matrix for the LLL 

Figure [S] shows numerical results for various matrix 
elements g a p(y) inside the LLL (0 < v < 1). The cal- 
culation is for the disorder potential shown in Fig. |3 
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for a magnetic field B = 3T. Two facts are immedi- 
ately apparent. First, for low filling factors [y < 0.34 
for this disorder realization), the conductance matrix is 
symmetric to high accuracy, g a p — gp a . From Eq. 
it follows that pp^ a — i-e. the electron has the 

same (very small) probability to scatter in both direc- 
tions. This is expected, since we know that at low filling 
factors, the states in the LLL are localized. Scattering 
from one terminal to another is only possible through 
tunneling of the electron from a state localized near one 
terminal, into a state localized near the other terminal. 
Such tunneling amplitudes are very small, explaining the 
small total probabilities for such scattering. Tunneling is 
also time-reversal symmetric - it happens with the same 
probability in both directions, consistent with the sym- 
metry of the conductance matrix. These ideas will be 
made more precise in Section where we analyze the 
general structure allowed for the conductance matrix. 

The second observation is that at high-filling factors 
v — > 1, we have g a , a +i — > e 2 /h (where if a = 6, a + 1 = 
1) while all other off-diagonal matrix elements become 
vanishingly small. In other words, an electron injected 
through terminal a + 1 scatters with unit probability into 
terminal a. This behavior clearly signals the appearance 
of the edge states at higher filling factors. These are chiral 
states, carrying electrons only in the direction consistent 
with the orientation (sign of) the magnetic field B. These 
states are localized on the edges of the sample, on the 
"vertical walls" created by the confining potential. 

Using Eq. (JSJ, this shows that in the limit v — > 1, the 
conductance matrix converges to the simple form: 



-1 1 



m ^ 5 (0) 



(-1 1 
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-1 1 
-1 1 









-1 1 
-1/ 



(18) 



The resistances corresponding to this asymptotic limit 
are straightforward to find. Solving /14 = g^ ■ V for 
I = 1,V 4 = 0, we find V 5 = U 6 = h/e 2 , V 2 = V 3 = 
0. Eqs. (133), (HQ) then give R" ifi2 = ^14,53 = V e2 ; 
^14,23 = -^14,65 = 0- I n other words, the values expected 
for the first IQHE plateau, which indeed is established 
when the LLL is more than half-filled. 



B. Fluctuations of resistances in the LLL 

The fluctuations of the various resistances are due to 
variations of g(v) from the asymptotic limit g^ . Using 
g(v) shown in Fig. |5j we calculate the resistances for all 
filling factors < v < 1. In Fig. we plot the pair i?^ 65 
and R\ A 53 as a function of v near half-filling. Three dif- 
ferent regimes are apparent: for v > 0.46, R H = h/e 2 
and R L — 0, corresponding to the first IQHE plateau. 
For 0.42 < v < 0.46, R L exhibits large fluctuations, 




0.2 0.4 0.6 0.2 0.4 0.6 



FIG. 5: Representative conductance matrix elements, in units 
of e 2 /h, as a function of v. The left panels show t/23, <?62, 345, 
(/i4, while the right panels show (732, 326, 354, 341- The top 
three left (right) panels characterize transport in (against) 
the direction of the edge currents. The bottom panels show 
the negligible direct transport between the current source and 
drain. For each pair, the traces are almost identical on the 
left, but very different on the right of the vertical line. 



however R H is still well quantized. This is precisely the 
type of behavior observed experimentally^ (see Section 
ITT1 point B). For v < 0.42 the transition to the insulating 
phase occurs, and all resistances increase sharply as the 
wave-functions at the bottom of the LLL become more 
and more localized, leading to progressively smaller val- 
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0.35 0.4 0.45 



v 

FIG. 6: A pair of R L and R H for the disorder shown in Fig. 
0and < v < 1. R L shows strong fluctuation before quanti- 
zation of R H is destroyed. The data are calculated at T = 0. 

ues for the matrix elements of g (see Fig. |SJ| . This results 
in large fluctuations of both resistances at T = 0; these 
are smeared out at finite T, as shown in Fig. [7| 

Figure shows all four resistances, for a different disor- 
der realization, at T = 11.6mK. The two Hall resistances 
are quantized to h/e 2 down to v « 0.4. Below v « 0.46, 
both R L exhibit fluctuations which quickly evolve into 
high peaks around v = 0.42 where both R H still have 
only minor deviations from h/e 2 . Comparing the upper 
and lower panel of Fig. [7| we see that the fluctuations 
of i?f 4 23 and i?f 4 65 are almost identical in this regime. 
Below v « 0.4 is the transition to the insulating phase. 

The resistances in the insulating phase in Figs. [S] and [7| 
are not perfectly similar to the experimental results^* In 
the experiment, even at the lowest temperature, the rise 
in R L is less steep and less noisy than the behavior shown 
by our simulations. This is because when the Fermi en- 
ergy resides amongst the low-^, strongly localized states 
of the LLL, the charge transport described here is no 
longer the dominant mechanism for conduction. Tem- 
perature assisted hopping conductivity and other inelas- 
tic scattering processes are likely to contribute to the 
conductance matrix more than the "band" conductance 
matrix which the simulation actually calculates. There- 
fore, we do not claim to have a realistic description of 
the insulating phase here. However, closer to half-filling 
the matrix elements of g become reasonably large, elas- 
tic scattering becomes the dominant mechanics for con- 
duction and our simulations describe well the behavior 
observed experimentally. 



C. The conductance matrix of higher LL 

The conductance matrix when the Fermi level is in a 
higher LL can be obtained from the conductance matrix 
of the LLL using the superposition principle,-- We simply 




0.48 



V 



FIG. 7: Both pairs of R L and R H for a different disorder. 
Both R L show strong fluctuations before the quantization of 
the R H disappears. The results correspond to T = 11.6mK. 

add a contribution g^ ' for each completely filled LL to 
the matrix g{v) describing the response of the LL hosting 
the Fermi energy (in this context, < v < 1 denotes the 
filling factor of the LL hosting the Fermi energy, not the 
total filling factor). As discussed, g^ describes transport 
through the edge states, which is the only possible con- 
tribution of a completely filled LL. Indeed, we see that if 
Ep is such that n LL arc completely filled, according to 
the superposition principle we have g = ng^ . This im- 
mediately leads to the solutions R H = h/{ne 2 ), R L = 0, 
i.e. the n th IQHE plateau, as expected. 

The validity of the superposition principle (which we 
take for granted here), is based on the belief that the 
IQHE transition can be described within one single LL, 
and that each plateau-to-plateau transition is in the same 
universality class. Wei et a/i experimentally confirmed 
this hypothesis. Shahar et al*Q> further developed this 
idea by mapping the transition between adjacent IQHE 
plateaus to the insulator-to-QH transition in the LLL. 

The superposition principle relies on several condi- 
tions, which are satisfied in our simulations: (1) One 
must confirm that for a single LL, as the Fermi energy 
is raised and v — > 1, the conductance matrix g — > g^ at 
reasonable filling factors, close to the center of the LL. 
We have already shown that our simulations satisfy this 
for the LLL. Higher LL should behave similarly, since the 
only difference is a simple change in the matrix elements 
of the disorder plus confining potential [see Eq. i|15|) and 
discussion following it]. (2) The cyclotron gap is large 
enough to justify neglect of the LL mixing. This is cer- 
tainly valid for v ~ n + 0.5, because here the charge 
transport is dominated by processes in the bulk of the 
Hall bar, where the disorder is much smaller than the cy- 
clotron gap. However, for the spatially close edge-states 
of different LLs, one does expect inter-Landau level and 
spin-dependent scattering to become relevant^ As a re- 
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FIG. 8: A typical transition from the first to the second IQHE 
plateau at T = 0. Lower panel: Ri4^2 an d ^14,23 in units of 
h/e 2 . Upper panel: The sum Ri4 :62 + R?i,23 = ^14,63 of 
the resistances shown in the lower panel, and R 2 t = R&3.63 
(shifted by — 0.5/i/e 2 ). Vertical lines indicate the boundaries 
of the critical region. See text for more details. 



suit, we do not claim that the edge-states produced in 
the single LL simulation are necessarily realistic. How- 
ever, this has no influence on the conductance matrix. 
The reason is simple: there is no backscattering among 
the edge-states £ All edge-states near the same boundary 
of the sample carry currents in the same direction, ir- 
respective of their LL index and spin-polarization. Scat- 
tering among these edge-states results in a unitary trans- 
formation which conserves the total probability to carry 
electrons from one terminal to the next. As a result, 
even though inclusion of the filled LL states in the sim- 
ulation may change some details of the edge-state wave- 
functions, the total conduction matrix g is not affected. 
This, in conjunction with the significant economy of CPU 
time when LL mixing is neglected, explains why we per- 
form the simulations for a single LL. 



D. Fluctuations of resistances in higher LL 

Using the superposition principle, the fluctuations of 
the resistances for 1 < v < 2 can be calculated using 
the conductance matrix g(v — 1) + g^°\ where g{v) is 
the LLL conductance matrix already analyzed and the 
second term is the contribution of the filled LL. 

The resistances R 4i 62 and i?f 4 23 corresponding to the 
disorder shown in Fig. [3] are plotted in the lower panel 
of Fig. [S] The three regimes found experimentally (see 
Sec. m points Al - A3) are clearly observed within the 
transition from the first to the second IQHE plateau. 
Their (approximate) boundaries are marked by vertical 
lines in Fig. |SJ as a guide to the eye. At low-^ (high-B), 
the fluctuations of R H and R L are correlated such that 
R L + R H = h/e 2 . This is seen more clearly in the upper 
panel, where their sum R^ 4 62 +-Ri4 23 = ^14,63 is plotted. 
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FIG. 9: The result of a simulation performed at T — 1 1.6mK. 
(a) Comparison between R 2 t and R H + R L . Here, Rit has 
been shifted up for clarity, (b), (c) the two pairs of R L and 
R H for 1 < v < 2. The vertical lines indicate the boundary 
of three distinctive fluctuation regimes. 



At high-!/ (low-S) R H = h/(2e 2 ) is quantized while R L 
exhibits strong fluctuations. In the intermediate regime, 
both R H and R L have strong, uncorrelated fluctuations. 
The upper panel of Fig. |H1 compares R L + R H with i?2t 
(i?2t is shifted by 0.5/i/e 2 ). As found experimentally, 15 
(Sec. [H] point C) the two curves are almost identical. 
Only minor differences at high-^ are visible. 

FigureElshows a set of data calculated for T — 11.6 mK 
with a different disorder realization. Both pairs of resis- 
tances are shown (the lower panels) . The thermal average 
was performed using Eq. JjJJ . At finite T, the traces are 
much smoother than at T — 0. The scale of the fluctu- 
ations in Fig. [5] is comparable to that of the experimen- 
tal data in Ref. - an indication that our simulations 
provide a good approximation of the experimental condi- 
tions. While the fluctuation patterns are clearly sample 
(disorder) specific, the three distinct fluctuation regimes 
are again clearly seen. The comparison between i?2t and 
-Ri4,63 (upper panel) shows that they are approximately 
equal over the entire transition. Moreover, it is also ap- 
parent that the two pairs of R H and R L have different 
fluctuation patterns, except on the high-z/ side where the 
two R L are almost identical. Another observation in all 
our simulations is that everywhere during this 1 < v < 2 
transition, < R L < 0.5h/e 2 and 0.5/i/e 2 < R H < h/e 2 . 
Within experimental error bars, the experimental resis- 
tance traces also stay within these limitsii 5 " Thus, all ob- 
servations are in agreement with the experimental facts. 
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V. GENERALIZED LANDAUER-BUTTIKER 
MODEL 

So far, we have shown that our numerical simulations 
faithfully reproduce the experimental results. In this sec- 
tion we explain the physics underlying these various types 
of correlations, using a simple but very general model. 
The consequences of changing the orientation of the mag- 
netic field, B — > —B, are also addressed. 

Since all resistances are functions of the various con- 
ductance matrix elements, it is clear that the correlations 
of the various resistances are encoded in the structure of 
the conductance matrix g. Thus, it is important to un- 
derstand the allowed structure of the conduction matrix, 
based on the various constraints it must satisfy, and the 
physics known to be relevant for transport in a LL. 

As already discussed, the conductance matrix in the 
presence of n filled LLs is of the form ng^ + g(y — n), 
where the first part is the edge-state contribution of the 
filled LL levels, and the second term is the contribu- 
tion of the LL hosting the Fermi energy. The challenge, 
then, is to understand g(v). Its matrix elements satisfy 
the following general constraints: (i) since they are pro- 
portional to transmission probabilities, all off-diagonal 
matrix elements are positive numbers g aj p^ a > [see 
Eq. (JZJ)]; (ii) Let Nt > 2 be the number of terminals 
(Nt — 6 is the case of interest to us). The Nt(Nt — 1) 
off-diagonal matrix elements must satisfy the 2Nt con- 
straints of Eq. JSJ. Half of them fix the value of the 
diagonal matrix elements, leaving a total of Nt — 1 sup- 
plementary constraints for the off-diagonal elements (one 
constraint is trivially satisfied if the other 2Nt — 1 hold). 

If Nt = 2, this implies immediately that the conduc- 
tance matrix must be symmetric, i.e. the most general 
possible form is: 



92 = G 



-1 1 
1 -1 



(19) 



G is the total conductance between the two terminals, 
since in this geometry we must have I\ = —1% = —I 
(all the current injected through one terminal must be 
removed through the other terminal), and therefore 



92 



Vi 
V 2 



I = G(Vi - V 2 ) 



For Nt — 3, the situation is more interesting. First, we 
separate the "symmetric" contributions. Let us denote 

G a (3 = min. (gap, 9pa) 

We can then rewrite: 

1 1 \ / -1 1 

.93 = G u ( 1 -1 + Gi 3 
0/ \ 1 -1 



+ Ag 3 




If the system is time-reversal symmetric, then A53 = 
since transport must proceed with equal probability in 
any two opposite directions, p a ^p = Pp^a- However, in 
the presence of a magnetic field, there is a preferred di- 
rection of charge transport chosen by the magnetic field. 
Thus, in cases of interest to us the conductance matrix is 
generally not symmetric (as already exemplified by g^). 

The matrix A53 which contains the non-symmetric 
contributions must still satisfy the constraints of Eq. JSJl, 
since the symmetric contributions do. Moreover, by con- 
struction, exactly 3 of its 6 off-diagonal matrix elements 
are zero. The remaining (positive) matrix elements must 
be all equal, and such that there is one on each row and 
on each column, otherwise the constraints cannot be sat- 
isfied. The only possibilities are either 



Ag 3 




or A33 




where c > is a constant. We call such contributions, 
which involve a closed loop of at least three terminals 
(here, 1^2^3orl^3^2)inan order selected by 
the magnetic field orientation, a "chiral" contribution. 

This approach can be straightforwardly generalized to 
Nt > 3. In cases with broken time- reversal symmetry, 
we expect to have some symmetric and some chiral con- 
tributions. To be more precise, we define the matrix 

l{ a , b)\ a p = SaaSpb ~ ^S aa Sj3 a ~ -^ab&Pb (20) 

which contributes a unit to the off-diagonal element g a b- 
For any ordered sequence of terminals a±, a 2 , a n , we 
define the matrix: 



r(ai, • • • 



,a n ) = l(a 1) a2) + i(a 2 ,a 3 )-\ \-l(a n ,ai). (21) 



Any such f matrix satisfies all the constraints of Eq. |JSJ. 
With this notation, and assuming that the magnetic field 
is such as to select the first form of we can rewrite: 



.93 = G 12 f(l, 2)+G 13 f(l, 3)+G 23 f (2, 3)+cr(l, 2, 3) (22) 

Any Nt x Nt conductance matrix g can be decom- 
posed in a similar fashion, into a sum of symmetric ( "re- 
sistances") and chiral contributions: 



9 = z2 c °i>- < a ^( ai > " ' ' a ™)- 



(23) 



where c ai) ... jan are positive numbers. The symmet- 
ric part is simply ^ Q<6 G a bf(a, b), where G ab = 
mm(g a b, gba)- This leaves at most Nt(Nt — l)/2 fi- 
nite, positive off-diagonal matrix elements, which can be 
grouped into a sum of "chiral" contributions. To do this, 
start with the smallest non-zero matrix element left, say 
gab — G a b, then draw a loop a — > a\ — > ■ ■ ■ — > a n — > b with 
bonds only connecting pairs of contacts sharing non-zero 
matrix elements. We can now separate a contribution 
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FIG. 10: A semi-classical illustration of possible transport 
between terminals, for (a) low filling factors v — > 0. In this 
case, all states in the LL hosting the Fermi level are local- 
ized, and transport can only occur by tunneling, (b) at high 
filling factors v — ► 1, the edge states are established and domi- 
nate the transport. Jain-Kivelson tunneling through localized 
states inside the sample can create short chiral currents. 

{{Jab — G a b)f(a, ai, a n , 6), insuring that the number of 
zero off-diagonal matrix elements of the remaining ma- 
trix has increased by at least one, while all the finite 
off-diagonal matrix elements are still positive. The pro- 
cedure is repeated until the decomposition is completed. 



A. Six-terminal geometry for the IQHE 

While the decomposition of Eq. (|23[) is very general, 
which of the decomposition terms are important depends, 
of course, on the physical system considered. We now 
return to the problem of interest to us, namely the 6 x 
6 matrix g(v) describing the contribution to the total 
conductance of the LL hosting the Fermi energy. From 
now, < v < 1 stands for the filling factor of the LL 
hosting the Fermi energy. The total filling factor is n + v, 
where n is the number of filled LLs. 

At low-v, all states of the LL hosting Ep are localized 
and transport can only occur through tunneling. Tun- 
neling occurs with equal probability in both directions. 
Thus, on very general grounds, here we expect the con- 
ductance matrix to be symmetric, with small off-diagonal 
elements since tunneling probabilities are small. 

Let us consider this case in more detail. In Fig. HUf aL 
we sketch some possible routes for charge transport be- 
tween three terminals, at \ow-v. The wave-functions are 
drawn in a semi-classical manner, following equipoten- 
tials of the disorder potential in the direction determined 
by the sign of B. This does not mean that our arguments 
only hold in the semi-classical regime. They are general 
and hold in the quantum regime - we just do not know 
how to draw quantum mechanical wave- functions. 

Most of the electrons injected through any of the ter- 
minals will be scattered back. However, with a small 
probability P12, electrons injected from terminal 2 can 
tunnel to near terminal 1. Another possible route of 
scattering from 2 into 1, is for electrons to first tunnel 
to a state near 6 (probability P25) and from there to near 
1 (probability p 16 ). The total probability for this process 
is (1 — pi2)p26P6i(l — Pii)- The electron can, however, 



make any number of closed loops between finally entering 
in 1, and so the total probability to arrive from 2 to 1 is 

(1 -pi 2 )p26P6i(l -P12) 

P2~>i = P12 H : 

1 - P12P26P61 

Similarly, an electron injected through 1 can scatter into 
2 either directly [with probability (1 — Pie)pi2(l — P26)] 
or can make any number of loops, yielding 

(1 -Pie)pi2(l -P2e) 

Pl^2 = 

1 - P12P26P61 

The other probabilities pi^e,Pe^i,P2^6 an d P6->2 can 

be calculated similarly. 

2 

Using Gab = \ mm(p a ^ b ,p b ^ a ) > 0, we find G a p = 

\Pafi + 0{p i ), i.e. indeed, the largest contribution to 
transport between these 3 terminals comes from direct 
tunneling and is symmetric. However, the conductance 
matrix is not fully symmetric, since P2^i — Pi^2 > 0. 
In fact, one finds that p2-*i — Pi^2 = P&^2 — p 2 ->6 = 
pi^ 6 - p 6 ^i = c = P12P2& + P12P16 + P2ePei + 0(p 3 ), 
showing also the appearance of a very small chiral contri- 
bution. Thus, the contribution of these processes to the 
conductance matrix is precisely of the expected form of 
Eq. (1221 . This derivation completely ignores interference 
effects between different scattering paths, and therefore 
is valid only in the presence of significant dephasing. Our 
numerical simulations, on the other hand, assume full co- 
herence between all electron wave-functions (there is no 
inelastic scattering in our Hamiltonian) . In such a case, 
one should sum the amplitudes of probabilities for various 
processes, and then square its modulus to find the total 
probability. The derivation for this case is very similar to 
the above one (also see Ref . Il7j) . The main contribution 
to the off-diagonal conductance matrix elements are still 
due to the direct tunneling, e.g. Gyi — ^- 1 ii2 1 2 , where 
ti 2 is the amplitude of probability to tunnel from 1 to 2 
(so that I i 1 2 1 2 = P\2). The difference is that the chiral 
current is now of order |i| 3 , not p 2 — |i| 4 as when the 
interference is ignored. The reality is in between, since 
in the real samples there is some amount of decoherence. 
Irrespective of how much decoherence there is, a pre- 
dominantly symmetric conductance matrix is inevitable 
if tunneling is the only means of charge transport. 

Of course, in order to derive an expression for the en- 
tire 6x6 conductance matrix, we have to also consider 
tunneling to the other 3 terminals in all possible combi- 
nations. It should be apparent that as long as all p <C 1, 
the only effect of that is to add symmetric terms of order 
p between all pairs of terminals connected by direct tun- 
neling, and much smaller chiral terms, of order p 2 or |i| 3 , 
for various closed loops. This explains why for low-i^, g is 
symmetric with small off-diagonal components, as indeed 
confirmed by the numerical simulations (see Fig. [SJ ■ 

At high-i/, on the other hand, the transport mechanism 
is very different. As expected on general grounds (and 
confirmed by the numerics) edge states are established 
within the LL once v > 0.5. With very high probability, 
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electrons are transported through the edge states to the 
next terminal a + 1 — > a. In the limit v — ► 1 we expect 
and find g -> ^f(l, 2, 3, 4, 5, 6) = g (0 \ 

For intermediate v, however, shorter chiral loops con- 
taining edge states can be established through tunneling, 
as sketched in Fig. HOf b). Assume that an electron leav- 
ing contact 3 can tunnel with amplitudes of probability 
£3 and t§ to and out of a localized state inside the sam- 
ple, to join the opposite edge state and enter 5. This is 
precisely the Jain-Kivelson phenomenology^ Summing 
over all possible processes and ignoring decoherence, we 
find their result 22 



P3- 



h 



tat 



3' 5 



1 - r 3 r 5 exp(i 



(24) 



where (j> is the flux enclosed by the localized state, 
0o = h/e is the unit of flux, and |r*3 j = yl — | ^ 3 1 2 , 1 7*5 j = 
y/l — jtgp are the amplitudes to avoid the corresponding 



tunneling. On the other hand, ^5^3 = 0, since in this 
scenario, no electron leaving 5 enters 3. Thus, in this 
case there is no symmetric term proportional to f(3, 5), 
only a chiral loop term 353^(3, 4, 5). Physically, this term 
represents the backscattered current of the Jain-Kivelson 
models Other shorter chiral loops can be established by 
tunneling between other pairs of edge states, so in the 
high-i/ limit we expect the conductance matrix to be a 
sum of such chiral loops. Again, inclusion of partial or 
total dephasing does not change this conclusion. 

This analysis has reconfirmed our assertion that the 
conductance matrix can be decomposed into a sum of 
symmetric terms and chiral terms. For this particular 
system, we have now shown that at low-i/, the symmet- 
ric terms are the dominant (although small) contribution, 
while at high-z^, chiral loops are the dominant contribu- 
tion. Near half-filling, we expect both types of terms to 
be important. Consider then the general form: 



g = ngW + G 12 f (1, 2) + G 16 r(l, 6) + G 26 f (2, 6) + G 34 r(3, 4) + G 35 r(3, 5) + G 45 f (4, 5) + 

c 5 (0) + cif (1, 2, 6) + c 2 f(2, 3, 5, 6) + c 3 f(3, 4, 5) + c 4 f(l, 2, 3, 5, 6) + c 5 r(2, 3, 4, 5, 6) 



(25) 



r 



The first term is the contribution of the n completely 
filled lower LLs, while all other terms describe trans- 
port in the LL hosting Ep. Eq. I|25l) is not the most 
general possible decomposition - that would have 15 
symmetric and 10 chiral terms. In Eq. I|25|) we as- 
sume no tunneling (no symmetric terms) between the 
left and right sides of the sample. The largest such 
neglected terms are G 2 3r(2,3) and Gs6r(5,6). Numer- 
ically, we find G 2 3,G56 < 10~ 4 ^- [see Fig. where 
G 2 3 = min(<7 2 3, <73 2 ), etc.]. We analyze the influence of 
these neglected terms on the resistance fluctuations in 
the next section. The 4 chiral loops neglected all have 
one "diagonal" link between 2 and 5, or between 3 and 
6. As explained, chiral loops are important at larger fill- 
ing factors, because of Jain-Kivelson scattering between 
opposite edge states. Such scattering does not mediate 
direct transport between 2 «-> 5 or 3 <-> 6 (see discussion 
below), so these terms can be safely neglected. 

The model of Eq. (|25|l is thus a parameterization of 
the 6x6 matrix g with 12 independent parameters. At 
low is, the chiral terms Cj — > and the decomposition 
is dominated by the symmetric terms. At high v, the 
symmetric terms vanish Gap — > while the chiral terms 
are important, in particular cq — ► 1. Near half-filling, all 
terms are likely to contribute. Eq. I|25l) thus covers all 
possible filling factors < v < 1. Figure ITT1 sketches the 
elements contained in this model. The closed directed 
loops indicate chiral terms while the resistors indicate 
the existence of symmetric terms. Resistors suggest dis- 



sipation; as we show shortly, they play a role only during 
transitions between QHE plateaus. Close to integer fill- 
ing, where the Hall conductances are quantized and the 
sample is disipationless, g — ► ng^ and all symmetric 
("resistors") terms vanish, as do the shorter chiral loops. 
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FIG. 11: Terms contained in the general decomposition of 
Eq. 12511 . Symmetric terms are represented as resistors, while 
chiral terms are represented as closed directed loops. 
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B. The correlations of the resistance fluctuations 

With Eq. (J23), the equations Iu= g-V and I 63 — g-V' 
can be solved analytically and the various Hall, longitu- 
dinal and two-terminal resistances can be calculated in 
terms of these 12 parameters. The complete solutions are 
very long and we do not list them here. 

The following identity is found to hold: 



R14. 



Re 



1 



e z n + c + c 2 + C4 + c 5 



(26) 



rt 14,62 



^14.23 



Since i?63,63 = R2t, whereas i?i4,63 
^14,65 + ^14,53: tlh s means that R 2 t = R H + R L irrespec- 
tive of the value of the 12 parameters. In other words, 
this identity is obeyed for all v, which explains why it 
is observed in both experiment and simulation. (This 
identity is not observed in experiment for low magnetic 
fields. This can be ascribed to deviations from the IQHE 
regime, where some of the approximations we make here 
- ignoring the LL mixing, for example - fail.) 

In Eq. (|2l)|) . n + c + c 2 + c 4 + c 5 is the total chiral 
current flowing along the 6 — > 5 and 3 — > 2 edges. As 
discussed, at low-^ the chiral currents in the LL hosting 
Ep are vanishingly small: Cj = 0, i = 1, 6 (there are no 
edge states established yet, and tunneling contributions 
to chiral currents are of order t 3 or less, as showed in the 
previous section. Below v c , all t < 10 -2 , see Fig. |SJ). It 
follows that at low-f, R L + R H = h/(ne 2 ), explaining 
the perfect correlations in the fluctuation patterns of the 
two resistances, seen experimentally and numerically. 

In the high-;/ regime the transport is dominated by 
the chiral currents created by tunneling between opposite 
edge states, through localized states inside the sampled 
The most general possible situation is sketched in Fig. 1121 
pi, p 2 , and P3 are the total probabilities for tunneling 
between the pairs of edges, summing over all the pos- 
sible tunneling processes through all available localized 
states in the sample. Each such contribution fluctu- 
ates strongly as B (or v) is changed, since the magnetic 
flux enclosed by various localized states changes signifi- 
cantly [see Eq. I|24|l for the simplest possible expression 
for such a probability]. The various conductance matrix 
elements can be simply read off this figure; for example, 
J?5i2 = P2->i — 1 — Pi- After adding the contribution 
of the n filled LL, the total conductance matrix can be 
written as a sum of chiral terms, consistent with Eq. (I25|l : 



g = (1 - pi - p 2 - P3 )g (0) + - P 2 [r(l, 2, 6) + r(3, 4, 5)] 

+^ [p 3 r(l, 2, 3, 5, 6) + Pl f(2, 3, 4, 5, 6)] + ng^ . (27) 
h 

In this case, the equation I 14 = g ■ V is trivial to solve. 
We find: 



^14,62-^14,53- (n + 1)e 2> 



• 



FIG. 12: Possible Jain-Kivelson processes in the sample for 
high-f region. 



i.e. the Hall resistances are precisely quantized, whereas 



R 



14,23 



R 



1>2 



14,65 



(n + l)e 2 n + 1 — p 2 



Note that the result only depends on p 2 , i.e. on tunnel- 
ing between edge states in the central part of the device, 
between the 4 voltage terminals. This is further justifi- 
cation that a short sample is sufficient for the numerical 
simulation. Since p 2 has a strong resonant dependence on 
Ep (or is), it follows that the two R L fluctuate strongly, 
but with identical patterns. This is precisely what is ob- 
served in experiment and numerics, on the high-i/ side of 
the transition (see Figs. [7| and [5J and supports our as- 
sertion that the fluctuations in this regime are caused by 
Jain-Kivelson tunneling. In particular, if n = (transi- 
tion inside the spin-up LLL), R L can be arbitrarily large 
when p 2 — > 1 (see Fig. [7J . In higher LLs, the amplitude 
of fluctuations of R L is h/[n(n + l)]e 2 or less, as observed 
both experimentally and in our simulations (see Fig. |5J. 

These results are very interesting because they show 
that one does not actually need to know the conductance 
matrix in order to understand what correlations might 
exist between various resistances. All that is needed is 
to have some idea of its general structure, which can be 
inferred based on physical arguments. At low v the con- 
ductance matrix is the sum between ng^ and a symmet- 
ric matrix describing tunneling between various terminals 
(ignoring tunneling between the left and right sides, for 
the time being). The solutions of I14 — g- V for such 
a matrix always satisfy R L + R H = h/(ne 2 ), no matter 
what are the off-diagonal values of the symmetric com- 
ponent. It follows that this identity must be obeyed in 
experiments as long as the conductance matrix has this 
form, i.e. for filling factors low-enough that tunneling 
is the dominant transport mechanism in the LL hosting 
Ep. On the other hand, for high filling factors, physical 
arguments based on the appearance of the edge states 
and the Jain-Kivelson phenomenology lead to the con- 
clusion that the general form of the conductance matrix 
is like in Eq. I|27l) . In this case, we find that irrespective 
of the values of the parameters pi,p 2 and p 3 , the R H are 
quantized while the two R L fluctuate with identical pat- 
terns. Finally, to cover all possible v, we must include 
in g both types of allowed symmetric and chiral terms. 
In this case, we find that the identity R 2t = R H + R L 



14 



holds, irrespective of the values of the various parameters. 
Such arguments can be straightforwardly generalized to 
geometries with any number of terminals, allowing one 
to easily test in what cases are correlations expected on 
such general grounds. The behavior when the sign of B is 
changed can also be understood easily, as we show now. 



negligible compared to the symmetric terms. To first 
order perturbation in G23 and G56, we find: 



^14,63 

n 



1 [n(G 26 + G12) + F^G 



23 



n 2 [nF 2 + F 1 ] 
HG35 + G 45 ) + F 3 ]G 56 
n 2 [nF 4 + F 3 ] 



(28) 



C. Changing the orientation of the magnetic field 



where 



If B changes sign, the Onsager reciprocal relation 18 
reads 

g(-B) = [g(B)f. 

That this must be so is obvious for our generic conduc- 
tance matrix of Eq. i|25|) : the time-reversal symmetric 
tunneling is not affected by this sign change, but the flow 
of the chiral currents is reversed (equivalent with taking 
the transpose). The solutions of j 14 = g(—B) ■ v are then 
related to the solutions of In = g(B) ■ V by v 2 = Vq, 
V3 = V5, V5 = V 3 and vq = V 2 , provided that the same 
index exchanges 2 <-> 6, 3 <-> 5, are performed for all G a b- 
Terms not invariant under this transformation are G12, 
Gi6, G43, and G45. In the experimental setup, these four 
terms must be very small, due to the long distance be- 
tween source and drain, and their nearby contacts. (In 
the simulation, these terms are sometimes not negligi- 
ble because the simulated sample is rather short). If 
we set these 4 terms to zero and keep only the largest 
symmetric terms G26 and G35 in Eq. I|25|l . we find that 
i?f 4 2 3(-B) = Ri 4 65 (—B) and vice versa. In other words, 
the fluctuation pattern of one R L mirrors that of the 
other R L when B — > —B, as observed experimentally - 16 : 
Small violations of this symmetry observed experimen- 
tally at low-z/, are likely due to the perturbative correc- 
tions from the very small, non-invariant tunneling contri- 
butions proportional to G12 — Gi6 and G43 — G45. This 
is confirmed in the next section. 



D. Small corrections to correlations and symmetry 

The most significant terms neglected in the general 
decomposition of Eq. (|2*5|) are G23r(2,3) and G56r(5,6), 
which are due to tunneling between the neighboring volt- 
age probes. All other tunneling terms neglected are 
smaller than these, because they are between contacts 
further apart. Remember that even these two terms are 
very small; numerically we found that they are of order 
10 ^- or less (see Fig. |5J). We now investigate whether 
the inclusion of these terms violates significantly the var- 
ious correlations established in their absence. For sim- 
plicity, in the rest of this section we set e 2 /h = 1, in 
other words we measure all conductances in e 2 /h units. 

We first investigate the effects of adding the G23 and 
G56 terms on the \ow-v correlations R H + R L = l/n. 
We set all Ci = 0, since at low-i^ the chiral terms are 



Fi — G16G26 + G12G16 + G12G26, F 2 — G16 + G26 + G12, 

F3 = G34G45 + G3 4 G35 + G35G45, F4 = G34 + G35 + G45 . 

In the given experimental geometry, we expect that 
Gi2,Gi6 <C G26 and G34,G45 <C G35, since terminals 
1 and 4 are very far from 2 and 3, respectively 5 and 
6, and thus the tunneling probabilities must be much 
smaller. Within this limit, the previous expression can 
be further simplified to: 

D 1 G23 + G56 

^•44,63 ~ 3 

n n° 

In other words, the identity R H + R L = i?i4,63 = l/n 
should indeed be obeyed with high accuracy, as long as 
the direct tunneling between 2 and 3, respectively 5 and 
6, as well as the chiral contributions, are indeed small. 

We now consider the effect of adding the G23 and G 56 
terms on the identity R L + R H = R 2t . After expanding 
to first order in these two quantities, we find: 

„ R _ G 45 G 56 Gi 2 G 23 , , 

-n-14,63 — -^63,63 — Z "j r ' ' ' • (^J) 

A\ A 2 

where we define Ct = Xh=o Ci an< ^ 

A 1 = (n + GT-ci + G 35 )(n + GT-ci-C3)(n+co-l-C3-|-c 5 ) 

M = (n + GT-c 3 + G26)(n- + GT-ci-c 3 )(n-|-co-l-ci+C4) 

Since G45, G56, G12, and G23 are all small numbers, we 
have Ai,A 2 ~ n 3 > 1, and therefore this correction is 
also small for all values of v. 

Finally, we investigate how the B-reversal symmetry 
is perturbed by a small imbalance between G34 and G45 , 
respectively G12 and G\q. As stated before, R H (—B) 
and R L (—B) are calculated with g T . For simplicity, here 
we take G23 = G56 = 0, and find: 

R^ t65 (B) - R^ 23 (-B) = G34 ~ G45 , (30a) 

RiaM B ) - Ri4.6 5 (-B) = Gl2 R Gl6 - (30b) 

where 

B l = G u G 45 + (n+co+c 3 +c 5 +G u +G45)(n+C T -c 1 +G 35 ) 
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#2 = Gi2Gi 6 +(n+c +c 1 +C4+Gi 6 +Gi2)(n+CT-C3+G26) 

The largest contributions to B\ and B 2 come from 
n, G35 and G 2 q and explain why these corrections to per- 
fect symmetry of i?f 4 . 65 (B) and R^ 4 . 23 (—B) are small. In 
fact, we can see that the larger the various chiral terms, 
the smaller these corrections should be. This agrees with 
the experimental data, where the largest violations of this 
symmetry are observed between the neighboring large 
fluctuation peaks in i?i4 ; 63, on the low-z^ side of the cen- 
tral regime. From Eq. (J2EJ), we know that those sharp 
peaks are caused by fluctuating chiral components of the 
conductance matrix, such as 02,04, and C5. From Eqs. 
(|3*U)l . we see that when these parameters are large (peak 
value in i?i4 ; 63), the difference between i?f 4 . 65 (_B) and 
Ri4- 23 (— B) is suppressed. When i?i4 ; 63 is closer to 1/n, 
however, the chiral terms are smaller and the corrections 
to the difference becomes more noticeable. 

In all these correction terms, powers of n appear in 
the denominator. Since n is the number of underlying 
filled LLs, these corrections should be smaller in higher 
LLs and therefore the symmetries and correlations should 
be easier to observe. They have been experimentally ob- 
served in higher LLsji& but they are not so clear as in the 
first plateau-to-plateau transition. There are two reasons 
for this. First, the sample quality is rather poor and only 
the first transition is clearly observed—*^ Secondly, a 
smaller magnetic field is needed to reach higher LLs. As 
the cyclotron energy is reduced, inter-Landau level mix- 
ing and other effects which we did not consider in these 
simulations may start to become important. 



VI. SUMMARY AND DISCUSSION 

To summarize, in this study we show that first- 
principles simulations of the IQHE in mesoscopic sam- 
ples, based on the multi-terminal Landauer formalism 
appropriate for non-interacting electrons, recapture all 
the non-trivial correlations and symmetries recently ob- 
served experimentally. Moreover, we explain how these 
correlations and symmetries are direct consequences of 
the general allowed structure of the conductance matrix. 

Similar to the experiments, we find that the IQHE 
transition in higher LLs is naturally divided into three 
regimes. On the \ow-v side of the transition, the LL 
hosting Ep is insulating. If tunneling between left and 
right sides is also small, the fluctuations of pairs of 
resistances are correlated with excellent accuracy, i.e. 
R H + R L = h/ne 2 . This condition is obeyed if the typical 
size of the wave-function (localization length) is less than 
the distance between contacts 2 and 3. When the local- 
ization length becomes comparable to this distance, edge 
states begin to be established and the correlation between 
R L and R H is lost. On the high-z^ side of the transition, 
the edge states are established and responsible for most 
of the charge transport. However, localized states inside 
the sample can help electrons tunnel between opposite 



edges, leading to back-scattering like the Jain-Kivelson 
model. In this case, our simulations show that the two 
R L fluctuate with identical patterns, while the R H are 
quantized. Tunneling between opposite edges is likely 
only if the typical size of the wave- functions is slightly 
shorter than the distance between opposite edges. It is 
then apparent that the central regime in Figs. [H] and 
corresponds to the so-called "critical region" , where the 
typical size of the electron wave-function is larger than 
the sample size (distance between contacts 2 and 3, at 
\cw-v, or between 2 and 6 at high-z/). In these meso- 
scopic samples, the voltage probes act as markers on a 
ruler, measuring the size of the wave-functions at the 
Fermi energy. To our knowledge, this is the first time 
when the boundaries of the critical region are pinpointed 
experimentally. This opens up exciting possibilities for 
experimentally testing the predictions of the localization 
theory of IQHE. 

Conductance fluctuations at the IQHE transitions have 
been studied before by several authors. Cho and Fisher— 
and Wang et czi— have focused on the two-terminal con- 
ductance, with numerical simulations based on Chalker 
and Coddington's network model— Ando has numeri- 
cally computed conductances for two- and four-terminal 
samples^ using the Buttiker-Landauer formalism, but 
in the Green's function formulation that requires the dis- 
cretization of the sample to a lattice model— Their re- 
sults are consistent with ours, showing fluctuations in the 
resistances which are random, sample-dependent, and of 
order of h/e 2 . Stfeda, Kucera and MacDonald 2 » were 
the first to predict the relations R H + R L = h/(ne 2 ) and 
i?2t = R H + R L , using an analytical analysis of the Lan- 
dauer formula for a "two-terminal" sample, but which 
has the upper and lower sides of each terminal held at dif- 
ferent voltages. In effect, their current source and drain 
(leads 1 and 4) are moved to infinity, and connected to 
the central sample through edge states only. This is a 
special case of our general model of Eq. I|25|l . Buttiker— 
presented a detailed analysis of four-terminal systems, 
similar to our analysis of the six-terminal conductance 
matrix. He predicted the symmetry relations between 
conductance measurements which exchange the role of 
the voltage and current leads. 

However, it is essential to emphasize that neither two, 
nor four-terminal samples (with terminals at well de- 
fined voltages) can simultaneously measure both a Hall 
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FIG. 13: Edge states in a four-terminal measurement. Leads 
1 and 3 are the current source and drain, while leads 2 and 4 
are the voltage leads. 
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and a longitudinal resistance. This is obvious for a two- 
terminal measurement, which can measure a single two- 
terminal conductance i?i2,i2- For a four-terminal geome- 
try (see Fig. 1131) . consider the case of n filled Landau lev- 
els, v = n. We know that each LL will transport charge 
through a set of edge states, such as those sketched in 
Fig. 1131 The conductance matrix in this case is: 
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Taking leads 1 and 3 as the current source and drain, 
I = (—7,0,7,0) and assuming that lead 3 is grounded, 
one can solve I = §4 ■ V for the other 3 voltages. One 
finds trivially that V\ = V4 = Ih/(ne 2 ) and V2 = V 3 = 0. 
Thus, even though at first sight one might assume that 
one can measure a Hall resistance R H — (V4 — V2)// and 
a longitudinal resistance R L = (Vi — V3)//, it turns out 
that in fact both measurements give resistances exhibit- 
ing plateaus at quantized values h/ (ne 2 ). In other words, 
both results are related to a Hall resistance (a longitu- 
dinal resistance should be zero when v — n. Note that 
at transitions between the plateaus, the two resistances 
need not be equal, since the partially filled LL will intro- 
duce other off-diagonal matrix elements. Using the type 
of analysis we introduced in Section [V] one could now 
easily study what types of correlations might be possible 
in such a geometry). If one now compares this to a 6- 
tcrminal case, it becomes obvious that the four-terminal 
resistances have no reason to be related to the R H mea- 
sured in the six-terminals by experimentalists. Thus, in 
order to simulate and understand the experimental re- 
sults, it is essential to use the six-terminal geometry. 

These considerations show that for mesoscopic sam- 
ples, a full specification of the experimental setup is ab- 
solutely necessary for any interpretation of the measured 
quantities. For example, the various correlations and 
symmetries that we studied in this paper only appear 
in the six-lead setup. Other geometries can be analyzed 
similarly. A significant result of our work is the proof 
that an understanding of various possible (robust) cor- 
relations can be obtained based on simple arguments re- 
garding the general allowed structure of the conductance 
matrix, without need for detailed numerical simulations. 

Furthermore, we have shown that the full DC response 
function (the conductance matrix) of a mesoscopic sam- 
ple is characterized by a large number of parameters (12, 
for our six-terminal geometry, when we ignore the small 
tunneling between the left and right-side of the Hall bar). 
This is to be contrasted with macroscopic samples which 
only require two parameters, the Hall and longitudinal 
conductivities a xy and o~ xx , to fully characterize their DC 
response. A mesoscopic sample obviously has more de- 
grees of freedom to display fluctuations of its resistances. 



APPENDIX: CHOOSING THE RIGHT 
CONFINING POTENTIAL 

The confining potential is needed to define the Hall bar 
out of the larger area spanned by each LL Hilbert sub- 
space. It is convenient to make a symmetric choice. We 
take the confining potential to be negative in the region 
\-L x /2,L x /2) x (-£ y /2,0) (which is thus the Hall bar) 
and positive in the remaining region [— L x /2, L x /2) x 
(0, Ly/2) (the inverted Hall bar) and such that V c (x 1 y) = 
— V c (x,y + Ly/2). On the y-edges of the Hall bar we 
then have V c (x,—L y /2) = V c (x, 0) = 0. Inside most of 
the Hall bar, the confining potential is equal to — V^ ap , 
where V gap > is large enough to confine the electrons 
to the Hall bar when disorder is added, for all Ep < 0. 

The question is how should V c {x, y) behave near the 
Hall bar edges, at x = ±L x /2 or y = —L y /2, 0. For 
instance, V c could vary smoothly or abruptly from —V gap 
to 0, as the y-edges are approached. On the x-edges, we 
could take V c = or — V gap - Two of the possible choices 
are shown in Fig. The shape of V c near the edges is 
critical for the quality of the simulation, as we show here. 

Since little is known about the real shape of the con- 
fining potential, we test the various choices by analyzing 
the IQHE transition from R H = h/e 2 to R H = h/(2e 2 ), 
in the absence of disorder. The computation is done as 
described in the text, but with Vd — 0. We add a to 
the calculated g(v), to account for the contribution of an 
underlying, filled LL. These simulations were carried on 
smaller samples, to save CPU time. The results should 
show: (1) a sharp transition of R H from h/e 2 to h/(2e 2 ); 
(2) the built-in particle-hole symmetry of the Hamilto- 
nian, i.e. symmetry of the results about Ep = (this is 
a consequence of choosing a symmetric confining poten- 
tial and symmetric contact states for the terminals); (3) 
the two Hall resistances are always equal, and so are the 





FIG. 14: Two of the tested confining potentials. Left: This 
potential has discontinuous changes at y — and y = ±L y /2, 
and is smoothly connected to zero on the x-edges. Right: An 
improved potential with open ends and smooth y-edges. 
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FIG. 15: Hall resistances calculated for the confining potential 
shown in Fig. 1141 (right'). #14,62 (grey) and #14,53 (black), (a) 
The overview of the entire energy interval scanned by Fermi 
energy Ef- Here, #14,53 has been shifted up by one unit, (b) 
Amplified view of the very small oscillations close to Ef = 0, 
AR H = R H — j^j. (c) Amplified view of the transition from 
R H = 4 to # H< = A. 




E F (eV) E F (eV) 

FIG. 16: Longitudinal resistances calculated for the confining 
potential shown in Fig. I14f right ), #14,62 (grey) and #14,53 
(black) . (a) The overview of the entire energy interval scanned 
by Fermi energy Ef- Here, #14,53 has been shifted up by one 
unit, (b) Amplified view of the very small variations close to 
Ef = 0. (c) Amplified view of the IQHE transition. 



two longitudinal resistances (in the absence of disorder, 
the potential has rectangular symmetry.) 

Figure 1141 shows two confining potentials which were 
used for testing purposes. Typical results of tests using 
the potential in the right panel of Fig. El ar e shown in 
Fig. El and El We found several features of the confin- 
ing potential that result in unrealistic, unphysical results. 
First, abrupt changes of the confining potential near the 
y-edges, such as shown in the left panel of Fig. El result 
in strong, fast oscillations in the resistances (similar to 
those displayed in Fig. I15f b). but with an amplitude of 
order h/e 2 ). Secondly, raising the confining potential to 
near zero at the x-edges, as shown in the same case, also 
leads to undesired consequences, for example half quan- 
tization R H = R L = h/(2e 2 ) over a considerable energy 
range. This is because leads 1 and 4 are only connected 
to the left-most and right-most LL states. In this case, 
the energies of these states are close to zero, so they ap- 
pear to be huge barriers blocking the electrons on the lead 
to enter the sample. Electrons can only tunnel through 
these barriers with very small probabilities, and no chiral 
currents can exit at the rc-edges unless Ep w 0. In other 
words, this is equivalent to having very bad contacts, a 
situation which is avoided in the experiments. 

The right panel of Fig. 1141 shows an improved confin- 
ing potential. The confining potential varies smoothy 
near the y-edges, and the cc-edges are no longer fixed at 
constant value. The y-edges are smooth because in a 
real sample, we do expect a smooth rise of the confining 
potential on the edges of the Hall bar. Its gradient is 
determined by how the 2DES is defined on the sample, 
as well as various screening effects. 

Figure El plots the Hall resistances calculated with V c 
shown in Fig. El (right). Panel (a) shows an overview of 



the data. Here V gap = lOmeV, so that Ep was varied over 
a range spanning the entire confining potential. Particle- 
hole symmetry is obvious. Panel (b) zooms in on the 
central region close to Ep = in panel (a) and reveals 
the minute scale of oscillations in R H . In this regime, 
the edge states have been fully established, and the total 
conductance matrix is almost exactly 2g(°\ The small 
oscillations also exhibit the built-in particle-hole symme- 
try. Panel (c) zooms in on the transition region. One can 
see that the transition from first to the second plateau 
occurs within a small energy interval. This is because in 
the absence of disorder, the edge states are established 
immediately once Ep > — Vg ap . In all panels, we see that 
the two Hall resistances are indeed identical. 

Figure El shows the longitudinal resistances for the 
same simulation. All the expected symmetries are again 
observed. Except for the very narrow interval of energies 
where the transition takes place, the longitudinal resis- 
tances are vanishingly small. This is expected, since for 
all energies except near ±V gap all transport is due to edge 
modes, which do not contribute to R L . The small devia- 
tions from zero near Ep = are magnified in panel (b). 
In panel (c), we show the peaks associated with the tran- 
sition. The symmetries observed in the results confirm 
the numerical accuracy of the simulations. We are now 
confident that in the absence of disorder, the simulation 
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FIG. 17: A sketch of the sample used by Peled et aZ. t 14 i 15 i 16 
drawn to scale. The region spanned by our LL Hilbert sub- 
space is comparable to the size of the shaded area. 
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shows a clear integer quantum Hall transition. Thus, we 
have confidence in attributing extra features in the full 
simulations to the effects of the disorder. 

As we just showed, the confining potential of 
Fig. [21 (right) gives physically reasonable results. How- 
ever, we make one more modification, to account for the 
short size of the samples we use in numerical simulations. 
Fig.ll7lshows the geometry of the real sample used in the 
experimentsii^ii^iiii The length of the real sample is ap- 
proximately 20 /im. Given our computational resources, 
we simulate the relatively short shaded area, so that af- 
ter adding the confining potential, our Hall bar has the 
same thickeness and distance between the 4 central volt- 
age terminals, but is only 4 /im long. Thus, our current 
source and drain (leads 1 and 4) are much closer to the 
voltage probes than in experiments. 



Because of this short distance, we frequently observe 
direct tunneling between leads 1 or 4 and their near- 
est neighbor terminals (2, 6 respectively 3, 5). In the 
simulation, this kind of tunneling leads to large sym- 
metric matrix elements in the conductance matrix, e.g. 
.912 ~ .921 ~ e 2 /h, and abnormal behavior of R L and R H . 
Experimentalists take great care to avoid direct tunnel- 
ing between leads, which constitutes a short-circuit. To 
avoid such direct tunneling in the simulation, we add tri- 
angular potential barriers in the corners of the Hall bar. 
With these, the edges are effectively prolonged, corre- 
sponding to an increased separation between the voltage 
probes and current source and drain. Fig. shows the 
final form of the confining potential used in our simula- 
tions with six terminals. 



* Current address: Computer Science and Mathematics De- 
vision, Oak Ridge National Laboratory, P.O. Box 2008, MS 
6164, Oak Ridge Tennessee, 37831, USA; Center for Sim- 
ulational Physics, University of Georgia, Athens Georgia, 
30602 USA. 
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